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^ '. ABSTRACT 

fT^ I Terrestrial timescales show instabilities due to the physical limitations of the atomic 

CO . clocks. Stricter environmental isolation and increased numbers of improved cesium 

clocks and cavity-tuned hydrogen masers have resulted in time scales more accurate 
by a factor of about five. The use of different clock ensembles results in measur- 
ij-j I able changes in some millisecond pulsar timing data. We investigate the possible 

On ' application of millisecond pulsars to define a precise long-term time standard and 

positional reference system in a nearly inertial reference frame. Although possible 
quantitative contribution of the two longest studied millisecond pulsars to terres- 
' I trial timescales appears minimal, they may prove useful as independent standards 

^ . in identifying error sources that are difficult to detect due to the finite lifetime 

' and common reference frame of terrestrial clocks. New millisecond pulsars, per- 

. haps some with even better timing properties, may be discovered as a result of the 

^ ' current global pulsar search efforts. 

•i-H . 

^ ' 1. Introduction 

■ - - ' The 1982 discovery of the 1.6 millisecond pulsar B1937-|-2lEl has provided an 

object with rotational stability comparable to the best atomic time standards over 
periods exceeding a few years.0 The microwave beams of millisecond pulsars, which 
are rigidly anchored to rotating neutron stars, can act as precision celestial clocks 
owing to the combined effects of their large rotational energies, 10^°~^^ erg, and 
low energy loss rates. Millisecond pulsars provide an opportunity for astronomical 
observations to define International Atomic Time (TAI) over long periods. 

Because each pulsar has a different spin and spin-down rate, timing measure- 
ments of millisecond pulsars cannot provide a more accurate timescale than current 
atomic timescales. Since the best foreseeable daily-averaged pulsar timing measure- 
ments will be accurate to 0.1 microsecond, two observations separated by 5 x 10^ 
seconds would be needed to reach the 2.5 x 10~^^s/s TAI monthly precisiorM (de- 
fined as the internal rms of the data set, ignoring systematic errors), although more 



o 



frequent pulsar observations would reduce this considerably. 

The existence of a number of millisecond pulsars distributed across the sky leads 
to the possibility of timing several pulsars against each other, with terrestrial clocks 
providing merely the endpoint times and a means of interpolating observations.Q&Q 
Such a program can smooth out variations in TAI, better determine the masses of 
the outer planets,Q place liffUts on perturbations produced by a primordial spectrum 
of gravitational radiationjOB constrain the nature of interstellar turbulence,li3 and 
provide an astrometric tie for the planetary and radio reference frames.li^ 

The search for new millisecond pulsars is a matter of intense activity and steady 
progress. More than 30 millisecond pulsars have been found in globular clusters, and 
over 25 have been found in the galactic disk, over a wide range of galactic latitudes 
and longitudes. Unlike slow pulsars, more than 75% of the millisecond pulsars have 
stellar companions; millisecond pulsars are often termed recycled pulsars because 
it is believed they were spun-up by their companions.0 It is estimated that over 
30,000 in the galaxy are beaming in our direction.lli Non-cluster millisecond pulsars 
are discovered at a rate of about 1 per 200-300 square degrees.0 Assuming a factor 
of two improvement in search efficiency and ignoring serendipity, perhaps up to 400 
field millisecond pulsars could be discovered within the next few decades. 

In this work we summarize the problems and status of terrestrial timescales and 
pulsar observations, and explore the contributions pulsars can make to terrestrial 
timescales. 



2. Terrestrial Clocks 

According to the International System (SI), the second is defined so that the 
frequency of the (F, mj)=(3,0) to (4,0) hyperfine transition of cesium 133 is exactlv 
9.192631770 GHz, in the absence of a magnetic field and located on the geoid.t3 
As recognized in the definition, the great majority of the timing data used to com- 
pute TAI are from cesium-based frequency standards.li3 Cesium atoms in the (3,0) 
state are created in an oven and isolated from cesium atoms in other states by 
their deflection in an inhomogeneous magnetic field. They are then allowed to 
travel through a magnetically shielded cavity, in which they are exposed to 9 GHz 
microwaves at each end. The probability that the 9 GHz microwaves will induce 
a transition to the f4,0) state has a Ramsey-pattern dependence upon frequency 
of the microwaves;cJ the number of ions which have made the transition to the 
(4,0) state is measured by the current created after deflection by an inhomogeneous 
fleld at the exit point of the cavity. The principal error sources are the magnetic 
fleld correction and the phase difference between the microwave flelds at the two 
interrogation regions.El Recently, Hewlett-Packard introduced their model 5071 ce- 
sium standards, with a signiflcantly improved electronics packageJll This model 
cesium clock has demonstrated superior performance both in the laboratory^ and 
in the fleld.ill The world's most accurate laboratory standards are maintained by 
the Physikalisch-Technische Bundesanstalt (PTB) in Germany and the National 
Institute of Standards and Technology (NIST), which achieve long-term stability 



through use of special beam-path and field switching techniques and the use of a 
longer cesium travel path.c3 

A second frequency standard used for high precision time keeping is the hy- 
drogen maser. Hydrogen masers are based upon the 21-cm transition of atomic 
hydrogen and can achieve much greater precision in the short-term.c3 However, 
their long-term accuracy is limited by variations in the dimensions of the confining 
cavity and interactions with the cavity walls-E^ Special Tefion coatings have been 
developed which minimize the wall-shift. Recently the Sigma Tau Standards Cor- 
poration has begun marketing masers whose cavity capacitance is stabilized and 
optimized through an attached varactor diode and a tuning circuit which measures 
the maser amplitude at 7.5 kHz on each side of the microwave resonance.ii 

The precision of the best types of frequency standards at the U.S. Naval Ob- 
servatory (USNO) as a function of the averaging time are shown in Figure 1, which 
was adapted from Breakiron.i2l Note that masers have significantly better perfor- 
mance on short timescales, but on the longest scale, 60 days, theirperformance is 
similar to the cesiums. Following the methodology of Barnes et al.f3 we summarize 
in Table 1 the error performance of the clocks as a sum of white and random walk 
frequency components. Since these clocks have not been in existence long enough 
to gather data on still longer timescales, it is likely that future analyses will identify 
other components to the variation of the cesium standards. A long-term analysis 
of maser data will be more difficult to perform since they undergo discontinuous 
changes in frequency and frequency drift (the time derivative of the frequency) on 
timescales of a few months, which are removed by comparison with other clocks. 

Within the next decade, we can expect further improvements in terrestrial stan- 
dards. The closest to implementation are optically pumped mercury stored-ion de- 
vices. Their advantage is based upon the high frequency of the observed transition 
(40.5 GHz) and coherence times which allow line-widths as low as 17 mHz. Al- 
though the accuracy of prototypes was limited by vacuum contamination,!! recent 
designs based upon a linear ion trap have achieved frequency stabilities better than 
10~^^s/s on timescales of a few houraHa (frequency is defined hereafter as the time 
derivative of the timing error in units of time, not phase). Somewhat further away 
are clocks based upon trapped yttrium,@ single trapped mercury ions,0 and ce- 
sium atomic fountains, which achieve long coherence times by first cooling atoms 
in "optical molasses" to 700 nanokelvins and then allowing them to rise out of the 
measurement region until gravity brings them back.0 

3. Global Time Transfer 

The Global Positioning System (GPS) has within the last few years become 
the chief instrument for comparing clocks of different institutions, as well as for 
determining position on the Earth. It consists of a constellation of 24 satellites, 
run by the U.S. Air Force Second Satellite Operations Squadron, each containing 
at least one cesium frequency standard and circling the Earth twice per sidereal 
day.E^ These satellites constantly broadcast their coded position and timing infor- 



mation, in both classified and unclassified channels, enabling anyone with the proper 
equipment to determine, by triangulation, the position and time at the receiver. De- 
pending upon the kind of GPS receiver, the time can be communicated through an 
output signal and/or compared digitally to the time input from a standard at the 
site. The individual timing characteristics of the satellites are monitored by the 
USNO, which forwards the information to the Air Force for use in determining the 
frequency and time offset corrections and in measuring orbital variations. Although 
the specified GPS time precision is 1 microsecond, since January 1993 GPS time has 
kept within 260 nanoseconds of the USNO Master Clock #2 (UTC(USNO)), and 
the rms difference over that period was 60 nanoseconds. Satellites also broadcast 
on-line predictions of GPS— UTC(USNO), which are specified to be less than 90 
nanoseconds, but are usually below 20 nanoseconds (which is close to the measure- 
ment error). Evaluations of GPS— UTC (USNO) are published as USNO Series 4, 
and are also available via modem, Internet, and the World Wide Web home page 
http : //tycho .usno .navy .mil . 

While the received GPS signal is distorted and delayed by unmodeled tropo- 
spheric, ionospheric, relativistic, and multi-path effects, the largest error contri- 
bution to unclassified use is the intentional distortion of up to 340 nanoseconds 
imposed by the military for security reasons and termed "selective availability" 
(SA). SA consists of both the insertion of a variable time delay and a misrepresen- 
tation of satellite orbital parameters, which are used to determine the propagation 
time for the signal to reach the observer. 

If one is merely interested in time difference between remote clocks, it is possible 
to eliminate the error due to the variable time delay through the "Common View" 
method. In Common View, observations are scheduled so that both sites observe 
the same satellite at the same time. While the difference neatly cancels satellite 
timing errors, it cannot eliminate the errors due to the incorrectly broadcast orbital 
parameters. It is possible to correct for this after the fact, using satellite orbital 
elements routinely disseminated by the International GPS Service for Geodynamics, 
which relies upon multiple observations for its analysis. Common View techniques 
routinely provide time-transfer precision of 5 to 15 nanoseconds.il0 

While it is possible, in principle, to correct any GPS measurement in real time 
for SA using Common View, such scheduling tends to result in fewer observations 
which are also conducted at lower elevation angles, so that tropospheric and iono- 
spheric modeling errors are more significant. At the USNO, unclassified GPS trans- 
fers are achieved using the "melting-pot" method, which forsakes Common View in 
order to obtain a straight forward average of as many satellites as possible. Even 
in the presence of SA, a precision of 20 nanoseconds is attainable if one is^ble to 
average at least seventy 13-minute satellite passes over a 48-hour period.Ea The 
advantage of this method is even greater for receivers which can observe more than 
one satellite at the same time. Another form of time-transfer consists of using 
available non-GPS satellites to send signals to and from the two remote sites. Aside 
from achieving sub-nanosecond time transfer precision, such studies often reveal 
systematic and perhaps seasonal differences with conventional GPS measurements. 



of order 10 nanoseconds.il 

4. Terrestrial Timescales and the BIPM 

The goal of timescale formulation is to approximate "true" coordinate time, 
from an average of differential clock timing measurements. In essence, the clocks 
serve as frequency standards to create a heterodyned tone of known or defined 
value. Individual timing events are realized either through summation of cycles or 
use of "zero-crossings" of the sinusoidal output voltage. In practice, timescales are 
usually generated in the frequency domain, in which a mean frequency is computed 
and used to define the time from an arbitrary offset. 

Internationally, the official responsibility for timescale generation lies with the 
International Bureau of Weights and Measures (BIPM) in Paris, France as first 
established by the Convention du Metre in 1875. The BIPM produces TAI by in- 
dependently evaluating the individual characteristics of approximately 200 clocks 
maintained by almost 50 institutions.il Evaluations are made at 10-day intervals, 
and reported monthly in their Circular T. The BIPM also evaluates timescales gen- 
erated by many national institutions, including GPS time. Some of these timescales 
are free-running, such as the USNO "A.l" series, which is currently comparable to 
TAI in accuracy. Others, such as UTC(NIST) and UTC(USNO), are steered so as 
to track UTC closely. The UTC(USNO), which is the national standard for legal 
and military matters, is available on-line and is in practice based upon the output 
of a single maser that is steered daily through the addition of small frequency offsets 
so as to approach the computed mean of all contributing USNO clocks, which itself 
is adjusted towards extrapolations of the BIPM-computed TAI-UTC(USNO). As 
a result of this steering TAI and UTC (USNO) have since 1991 kept within 200 
nanoseconds. 

Once all the time comparisons have been assembled, the problem of timescale 
generation reduces to calculating weights and offsets for each standard. The BIPM 
first generates a free-running timescale, the Echelle Atomique Libre (EAL). The 
EAL is computed by evaluating clock weights each 60 days. The weights are deter- 
mined from the variance of the frequency offsets over the past year, after allowing for 
any intentional rate adjustments. A maximum weight is chosen so that the timescale 
will not depend upon just a few of the clocks. While both the general and specific 
procedures are periodically under review, in 1994 about half of all clocks used had 
this maximum weight, while most of the rest were weighted less than 20% of the 
maximum. The BIPM does not subtract offsets from the clock frequencies, and this 
can lead to minor discontinuities as clocks are added, subtracted, or reweighted. 
The TAI timescale is created by gently steering the EAL so that its long-term be- 
havior conforms to the "primary" frequencv standards, which are currently the two 
"long-tube" primary standards of the PTbH and NIST-7.i3 This is a practical way 
for a real time system to take into account the fact that the primary standards have 
greater long-term accuracy than accounted for by the precision-based weighting sys- 
tem. Again, there is a minor problem in that no allowance is made for fact that 



the NIST-7 frequency is corrected by about 1.9 x 10~^'*s/s for the AC Stark effect 
due to the microwave background experienced by the cesium atoms,E3 whereas the 
PTB frequencies are not. 

At the USNO, the A.l timescale is currently determined using a weighting 
function which weights clocks by their type. Recent maser data are given a much 
higher weight than recent cesium data, but older maser data are progressivelvdown- 
weighted, with zero weight assigned to maser data more than 60 days old.l23 Also, 
data from the older model cesium standards are included at 0.65 the weight of model 
5071 cesium clocks. Unlike the BIPM, each clock is introduced to the mean after 
subtracting an initial frequency offset. This has the effect of reducing discontinuities, 
and the drawback of subtly over-weighting initial clocks is soon lost in the random- 
walk noise as long as no systematic errors are present. For historical reasons, the 
A.l is intentionally steered a constant 19 ns/day£^ 

Using data from 1993 and 1994, the BIPM estimates the lOO-dav error of the 
EAL as slightly less than the TAI, and both to be about 2.5 x 10~^^s/s;B the errors in 
the 1980's were above 10~^^s/s. The recent improvement in both TAI and the USNO 
free-running timescale A.l is evident in Figures 2-4, which gives their difference as 
a function of time. Part of this apparent improvement is due to the fact that the 
contribution of the USNO clocks to TAI has varied (from 30% in the 1980's to a 
low near 15% in 1990 to almost 40% in 1994) but most is due to the world-wide 
increase in number, quality, and environmental isolation of the individual frequency 
standards. Similar increases in precision are shown in the comparisons with the 
independent PTB timescale. 

For practical reasons, other timescales are often used which differ from TAI by 
specific and calculable amounts. Coordinated universal time (UTC) differs from 
TAI only by an integral number of leap seconds, which crudely take into account 
variations in the rate of Earth rotation.Eil These leap seconds are extra seconds in- 
cluded in UTC at UT midnight on either December 31 or June 30 so as to minimize 
UTl— UTC, where UTl is tirne related to the rotational angle of the earth after 
allowing for polar motion.0' 00 They have lately been inserted at a rate approach- 
ing once a year, and on July 1, 1994, TAI— UTC= 29 seconds. Measurements and 
predictions of UTl, polar motion, and nutation parameters^ are available on-line 
from the home page pittp : / /maia . usno . navy . milj . In order to allow for continuity 
with ephemeris time (ET), Terrestrial Time (TT, also called Terrestrial Dynamical 
Time, TDT) is defined to be offset by exactly 32.184 seconds from the "ideal" form 
of TAI, which is what TAI would be if no instrumental noise were present. Although 
in the future TT may be realized from an average of many sources, including pul- 
sars, currently TT is realized only from TAI, by post-processing the xaw clock data 
and improving upon the initial weighting and steering correct ions. ^3 Since 1982, 
the BIPM-computed numerical difference TT— TAI has varied by over 3 microsec- 
onds absolutely, and over 2 microseconds after removing the quadratic component 
of the variation. In order to provide continuity in the interpretation of old data 
and allow for the effects of general relativity in a consistent manner, the lAU in 
1991 passed resolutions recognizing timescales which distinguish between centers of 



reference at the solar system barycenter (Barycentric Coordinate Time, TCB), the 
center of tha Earth (Geocentric Coordinate Time, TCG), and the geoid (TT, TAI, 
UTC, ET)BBi 

5. The Techniques of Pulsar Timing 

The precision of an individual pulsar time-of-arrival (TOA) measurement de- 
pends upon its pulse period and Jiie signal-to-noise ratio as determined from the 
particular instrumental setups.E^'H The fastest pulsar known today is still PSR 
B1937-I-21, whose 1.6 ms period is close to the theoretical limit determined from 
balancing the centripetal and gravitational forces, assuming a typical neutron star 
equation of state at densities of ^ 10^^ g cm~'^.BJ 

The simplest rotation model for a solitary neutron star can be given in terms 
of phase residuals as a power series in time: 

(j){t) =(j)o + 0(t - to) + <^{t - to)V2 + ^{t - tof/6 + ... (1) 

where (p is the pulse phase, O is the rotation frequency, O is the spin-down rate. 
The spin frequency is related to the pulse period by O = 27r/P (O = 27rz^), while 
O = —2nP/P'^. If the spin-down of the pulsar is driven by magnetic dipole radiation, 
due to misalignment of the magnetic and rotational axes, then the braking index 
n, defined by 00/0^ = 2 — PP/P^, is 3. If the initial pulsar period was Pq, the 
pulsar age can be computed from 

^ ^[l-(5)^"-^^]- (2) 



(n-l)P' P 

Most pulsars P measurements yield only upper limits, but in those pulsars for which 
P can be measured above the noise, n is found to be between 2 and 3. Inverting 
the equations above, we derive a theoretical value 

^ _ _ 2n7TP^ 

In the case of PSR B1937-I-21 the theoretical slowdown would result in a 0.1 mi- 
crosecond residual from a constant P model after 11 years. The relativistic time-of- 
flight effects and acceleration due to the galactic gravitational field can also make 
small, but non-negligible contributions to O and the theoretical Qxl Also, for 
pulsars in globular clusters, contributions to O from Doppler acceleration can be 
important where the gravitational potential of the cluster accelerates the pulsar, 
influencing the rotation period over several- year timescales.BEd 

Pulsar timing requires an explicit definition of both a steady atomic time stan- 
dard and accurate monitoring of the motion of the Earth with respect to the 
background stars. This task is accomplished by using a clock whose offset is fre- 
quently calibrated against a real time steered UTC service, such as UTC(NIST) or 
UTC(USNO). This real time UTC realization is later corrected to a more accurate 



timescale using after-the-fact determinations published by the BIPM. Although TAI 
is often used, TT is recommended as the most accurate terrestrial scale available 
to day .S Once a terrestrial timescale has been chosen, it must be converted to the 
barycentric frame (TCB). The omission of this time conversion would, among other 
things, result in a sinusoidal monthly error of about a microsecond as the Moon pulls 
the the Earth back and forth through the solar gravitational potential well. In most 
software, the time conversion happens automatically as part of the four-dimensional 
general-relativistic transformation which accounts for the Earth's orbit. The lAU 
has recommended that all users avoid transformations to Barycentric Dynamical 
Time (TDB),e3 which is scaled to remove secular differences with TT; as a result it 
diverges from TCB at a constant rate of 1.55 x 10~^, which would be absorbed into 
the fitted parameters. 

The geometric, general relativistic, and astrometric effects due to the position 
and velocity of the antenna with respect to the solar system barycenter must also be 
removed. In most cases this is done using the lERS values for UTl— UTC and the 
DE200 ephemeris.il Since most observations used for computing the Earth's orbit 
are Earth-based they are relatively insensitive to the uncertainties in the masses of 
the outer planets. The most uncertain long-term component of the ephemeris for 
pulsar work is the motion of the Sun and known planets about the solar system 
barycenter. The largest error contribution is due to the uncertainty in the mass of 
Saturnjwhich could result in a TOA error of order 7 microseconds with a 30-year 
period.E^I This and other such errors would be identifiable through their different 
and systematic projections upon the pulsars observed.Q Fukushima has examined 
the differences between the DE200, DE102, and DE245 ephemerides and found the 
differences between them (including their effects on the relativistic time correction 
TCB— TT) to be negligible for pulsar work.ca 

Pulsar signals must be corrected for dispersion caused by the intervening plasma 
using the frequency dependence of the delay. For a homogeneous and isotropic 
medium the group velocity is given by Vg = c(l — Wp/a;^)^/^, where tOp is the plasma 
frequency and co is the wave frequency. The plasma frequency in Gaussian units is 
Ljp = A-Knee"^ /m, where Ue is the mean electron density, and e and m are the electron 
charge and mass. The delay is proportional to the column density of electrons in the 
line of sight towards the pulsar. This column density is called the dispersion measure 
and is defined as DM = Jq Uedl, where z is the pulsar distance. The dispersion 
measure is a directly measurable quantity determined from the differential pulse 
arrival time between two frequencies. To first order, observing Vg gives a time delay 
between two frequencies as 

t2-h = ^(wo"' - oj7^)DM. (4) 
mc 

Observations of dispersion measures toward 706 radio pulsars give values that range 
from 1.8 to ~ 1000 cm~^ pcH For some lines of sight the turbulent electron 
population may cause different wavefronts to travel along slightly different paths, 
which can lead to delays that are not removable by the two-frequency formula. 



although some improvement can occur with multi-frequency observations. The 
unmo deled delay, which depends upon the total electron column density and the 
turbulent properties on the interstellar electron population, has been estimated in 
the worst cases to be ctLthe sub-microsecond level for radio frequency observations 
at 1.4 GHz and above.ty 

Another important consideration is that millisecond pulsars have steep radio 
spectra and are generally weak enough that observations at the largest telescopes 
are required in order to achieve an adequate signal-to-noise ratio. An estimate of 
the SNR achievable at current and future sites is provided in Table 2. 

Once the TOAs from a pulsar are assembled and iteratively corrected for known 
effects, several pulsar parameters are typically solved for. Uncertainty in the a priori 
pulsar period and spin-down rate are identified through their quadratic dependence. 
The pulsar parallax is identified from the TOA by its semiannual periodicity, posi- 
tion errors through their annual periodicity, and proper motion by a growing annual 
periodicity. Pulsar parameters associated with orbits about companions are more 
complex.0'il More exotic error sources, such as gravitational radiation and un- 
known planets, can be identified through comparison of delays with other members 
of the pulsar ensemble. If the TOA errors are modeled as a multi-pole expansion 
over the sky, timing errors are related to the monopole term, ephemeris errors result 
in a dipole dependence, and gravitational radiation would cause a quadrupole and 
higher order effect .a 

6. Stability of Pulsar Time Standards 

Using data on the millisecond pulsars B1937-f21 and B1855-(-09 made available 
by anonymous FTP from Princeton University the program TEMP0E3 was used 
to solve for the pulsar periods, spin-down rates, position, and proi)er motion. The 
residuals to the solutions are shown in Figure 5. As Kaspi et alB found with the 
same data, the residuals of PSR B1937-I-21 appear to display red noise (long-period 
correlations, in this case well-modeled as a non-zero P), whereas the residuals of 
PSR B1855-f09 do not. Since these two pulsars are close together in the sky, the 
red noise is considered to be intrinsic to the pulsar B1937-f21 and not due to the 
timescale or ephemeris errors. The timing noise for millisecond pulsars see 
follow the correlation with P found for slower pulsars by Arzoumanian et ai: 

As = 6.6 + 0.61ogP, (5) 

where the "stability parameter" Ag = log (^^|i>|(10^s)^^ with a scatter of ± one 
decade. Pulsar B1937-f21 has a Ag = —5.5. If this inverse correlation is valid for 
millisecond pulsars, then the PSR J1713+0747 would be a particularly promising 
candidate to have a very low Ag stability parameter ,11 with a red noise less than 
that of PSR B1937-I-21 by a factor of about 5. 

The effects of different input timescales are shown in Figure 6, which plots 
the difference between the residuals of parameter fits to PSR B1937-(-21, which 
differ only in the input timescale assumed. For these and subsequent solutions. 



the parameters were the pulsar period, period derivative, position, paraUax, and 
proper motion. In Figure 7 these difference solutions are double-differenced with the 
differences between the timescales, as made available by the BIPM, with a second- 
order term removed. The double- differences indicate how much timescale variation 
is soaked up into the fitted parameters, and hence is indeterminable from pulsar 
TOA data. The annual signature evident in the plots, for example, shows how pulsar 
timing data alone will never be able to identify strictly periodic annual variations 
in terrestrial standards. In the presence of an independent means of determining 
pulsar positions, proper motion, or parallax (such as astrometric Very Long Baseline 
Interferometry, or VLBI), such variations would be recognizable as displacements 
in the fitted parameters. Table 3 shows the current positional precision of a number 
of millisecond pulsars, derived from timing data. The coupling of the dynamic and 
radio reference frames will, along with VLBI observations of millisecond pulsars, 
aid in removing two degrees of freedom from the pulsar fitting models and improve 
our ability to constrain other sources of errors. 

Slightly extending the analysis of Blandford et a/. ,11 it is possible to apply 
a Weiner filter to the residuals by Fourier-transforming the residuals, and then 
transforming back to the time domain after multiplication by a "transfer function," 
which corrects for the distortion caused by fitting to pulsar parameters. The exact 
correction factor would also include a factor incorporating the clock and pulsar 
stabilities as a function of frequency. Since pulsar data are much less accurate on 
short time scales, compared to terrestrial standards, pulsar data would be given 
insignificant weights at short periods. For the generation of the plots in Figure 8, 
a simpler analysis was used in which all periods less than 400 days were discarded. 
Figure 9 shows the "transfer function" for two different time ranges, which indicate 
what frequencies in the raw data are absorbed by the fitted parameters .1521 The low 
values at periods of one year and six months indicate that any noise in the data with 
those frequencies will be absorbed into the fitted pulsar parameters. It is evident 
that the accumulation of another decade of data will bring about a considerable 
improvement in the ability of pulsar timing data to discern long-term variations. 

The contribution millisecond pulsar data can make to a timescale is entirely 
dependent upon the quantity and quality of the pulsar TOAs compared to the 
terrestrial standards. Included along with the simple clock model in Table 1 are 
parameters for the timing stability of PSR B1937-I-21, which is characterizable 
as a random walk in frequency, or perhaps an integrated random walk, which is 
even redder.!^ These numbers were chosen to be consistent with the residuals 
plotted in Figure 5, and they show that PSR B1937-I-21 is less accurate than a 
single model 5071 cesium standard. If the relation between timing stability and P 
is valid, then the pulsar J1713-I-0747 would be five times more accurate, and an 
ensemble consisting of several such pulsars could approach the accuracy of today's 
ensemble of terrestrial frequency standards. More important than the numerical 
evaluations is the fact that pulsars can provide an independent means of evaluating 
clock performance, and can serve to constrain the amount of extremely red noise in 
terrestrial clocks, which would be undetectable on the short timescales available so 



far. 
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Table 1: Simple Noise Model 





White 


White 


Random Walk 




Phase Noise 


Frequency Noise 


Frequency Noise 




nanosecond 


10-^5s/s(t/days)-°-^ 


10-i5s/s(t/days)+°-5 


cesium * 





25 


0.5 


maser * 





2 


1.0 


PSR B1937+21 ** 


100 





3 



* It is assumed that the times of the phase, rate, and frequency drift discontinuities 
have been identified, and their effects removed by least-squares methods. Such 
methods are required on monthly scales for masers, and much less frequently for 
cesiums. It has the effect of masking higher order variations. 

** Although today's best Arecibo data for PSR B1937+21 have formal TOA errors 
of 200 nanoseconds, higher precision will be achievable with the GBT and upgraded 
Arecibo. The random walk noise, computed_from Monte Carlo simulations, is nu- 
merically larger than modeled by Petit et a/.E3 



Table 2: Antenna Parameters at 1.4 GHz 





T 

sys 


Gain 


Bandwidth 


SNR 


rp :^ 
-'■ int 




K 


K/Jy 


MHz 


2 Polarizations 


Hours 


Arecibo, upgraded 


20 


11 


200 


150 


0.3 


Arecibo, current 


40 


8 


40 


24 


10 


GBT 


10 


1.9 


200 


50 


2 


Effelsberg 


20 


1.5 


200 


20 


14 


Nancay 


40 


1.5 


200 


10 


55 


Parkes 


20 


.64 


320 


11 


50 


Jodrell Bank 


20 


.64 


320 


11 


50 


NRAO 140' 


20 


.26 


120 


3 


800 



* Time required to observe PSR J1713+0747 with a TOA error of 100 nanoseconds 
assuming a flux of 3 mJy and a pulse width of 18 degrees. 



Table 3: Current Pulsar Astrometric and Timing Precision 



Pulsar 


Observation 


Timing 


RA 


DEC 


Name 


Duration 


Precision 


Precision 


Precision 




(years) 


(fxsec) 


(mas) 


(mas) 


PSR B1257+120 


2.6 


2.3 


0.4 


1.0 


PSR J1713+0747§ 


1.8 


0.4 


0.2 


0.3 


PSR B1855+09i 


6.9 


1.0 


0.07 


0.12 


PSR B1937+2li 


8.2 


0.2 


0.03 


0.06 


PSR J2O19+24250 


2.7 


3.0 


0.6 


0.9 


PSR J2322+2057izl 


2.3 


2.9 


1.0 


2.0 



Figure Captions 

Figure 1: Observed frequency stability of USNO standards as function of averaging time. Plotted 
are the average Allan deviation (rms of differences between adjacent averages, divided by root2) of 
the frequencies of 6 cavity-tuned masers and 12 model 5071 cesiums as a function of averaging time 
r. The error bars indicate the ensemble rms. 



Figure 2: Frequency difference between free-running timescales of TAI, PTB, and USNO (A.l) in 
fs/s. The USNO and PTB timescales are independent, but TAI is gently steered towards the PTB, 
whereas the USNO clocks' contribution to TAI has increased to 40%. Nevertheless, improvement 
in the last few years is obvious. The lowest plot shows the difference between TT and TAI, where 
TT represents TAI recomputed using hindsight corrections to clock weights, offsets, and drifts. A 
constant frequency offset has been removed from all plots. Data were obtain from the BIPM by 
anonymous ftp to address 145.238.2.2 The time range shown is from January 1982 to December 
1993. 



Figure 3: As in previous figure, except that the vertical axis shows the timing difference, with a 
quadratic term removed. Pulsars cannot contribute to the quadratic terms as their true periods and 
spin-down rates cannot be determined. 



Figure 4: As in previous figure, except the time range of data is shortened to include only the latest 
and best data, from 1990 to 1994. 



Figure 5: Timing residuals to pulsar parameter solutions for PSR B1937+21 and PSR B1855+09, 
from October 1984 to December 1992. 



Figure 6: Difference between timing residuals of PSR B1937+21. Solutions differ only in their 
reference timescale, and each plot is the difference between the solutions with the indicated timescale 
and with TT. The data are from October 1984 to December 1992 



Figure 7: Double-difference plots subtracting the TEMPO residual differences of previous figure 
from the terrestrial timescale differences. Since pulsars cannot determine variations below second 
order, a quadratic term was removed from the terrestrial differences. The data are from October 
1984 to December 1992. 



Figure 8: Upper part shows the residuals for a fit to pulsar parameters for PSR B1937+21. Lower 
part shows the same residuals after using a Fourier analysis to extract all periodicities shorter than 
400 days. No corrections have been made for spectral leakage or aliasing. The data are from October 
1984 to December 1992. 



Figure 9: Transfer functions for 165 observations equally spaced over 8 years, and for 330 observations 
equally spaced over 16 years as a function of period (1 /frequency). In this plot, a value of 1 indicates 
no noise is absorbed in the fit to pulsar parameters, and a value of indicates that all noise is 
absorbed. 
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